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Abstract This paper combines experimental data with 
simple mathematical models to investigate the influence 
of spray formulation type and leaf character (wettabil¬ 
ity) on shatter, bounce and adhesion of droplets im¬ 
pacting with cotton, rice and wheat leaves. Impaction 
criteria that allow for different angles of the leaf sur¬ 
face and the droplet impact trajectory are presented; 
their predictions are based on whether combinations 
of droplet size and velocity lie above or below bounce 
and shatter boundaries. In the experimental compo¬ 
nent, real leaves are used, with all their inherent nat¬ 
ural variability. Further, commercial agricultural spray 
nozzles are employed, resulting in a range of droplet 
characteristics. Given this natural variability, there is 
broad agreement between the data and predictions. As 
predicted, the shatter of droplets was found to increase 
as droplet size and velocity increased, and the surface 
became harder to wet. Bouncing of droplets occurred 
most frequently on hard to wet surfaces with high sur¬ 
face tension mixtures. On the other hand, a number of 
small droplets with low impact velocity were observed 
to bounce when predicted to lie well within the adhering 
regime. We believe this discrepancy between the predic- 
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tions and experimental data could be due to air layer 
effects that were not taken into account in the current 
bounce equations. Other discrepancies between experi¬ 
ment and theory are thought to be due to the current 
assumption of a dry impact surface, whereas, in prac¬ 
tice, the leaf surfaces became increasingly covered with 
fluid throughout the spray test runs. 


1 Introduction 


Pesticides are used for the management of weeds, in¬ 
sects or pathogens in most agricultural cropping sys¬ 
tems worldwide. They are generally applied as a spray 
with the aim of covering all or part of the target (path¬ 
ogens, insects, leaves or other plant parts). To be ef¬ 
fective at controlling the pest, the droplets must be 
retained on the target surface and not lost as spray 
drift, or deposited on the ground and other non-target 
sites within the sprayed area. Within this context, we 
present here experimental data on droplet impaction 
events using real leaves, commercial spray nozzles, and 
a variety of spray adjuvants. 

While our motivation is to study droplet impaction 
in the context of agrichemical spraying, another closely 
related application of this type of work is the spread¬ 
ing of pathogens from leaf to leaf by rain droplets after 
impaction. Splash from rainfall or overhead irrigation 
is important in the dispersal of particles from natu¬ 
ral surfaces, since secondary drops produced when rain 
droplets impact a surface can pick up pathogens and 
deliver them to another plant. Rain-splash is probably 
the second most important natural agent for disper¬ 


sal of spores after wind (Huber et al. 1997 Fitt et 


ah, 1989); splashing provides an effective short-range 


dispersal, particularly when the leaf canopy is dense. 
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whereas wind is the mechanism behind long-range dis¬ 


persal (Calonnec et al., 20131. In contrast to pesticide 


application, droplet sizes from rainfall are much larger 
in size. Typically the diameters of rain drops are greater 
than 1mm, whereas for pesticide application median 


diameters are typically less than 0.5mm (Dorr et al. 


20131. Further, in a rainfall event, leaves are typically 


wet and hence are more prone to shatter when com¬ 
pared to field application of pesticides, where most of 
the droplets impact upon a dry area of the leaf. Nev¬ 
ertheless, the shatter events which occur during either 
pesticide application or via rainfall are largely deter¬ 
mined by a competition between inertial and surface 


tension forces acting on the droplet at impact (Saint- 


Jean et al., 2006). 


There are three main possible outcomes when a drop¬ 
let impacts the target: it can adhere, bounce off the 
surface, or shatter into a number of smaller droplets 
(Mercer et al. 2010| . Each of these three outcomes can 
take different forms. If a droplet adheres it can return 
to the shape of the original droplet and rest on the sur¬ 
face, or it can spread over an area much larger than 


the initial droplet diameter (Xu et al. 20111. If a drop¬ 


let bounces it can fully rebound or part of the droplet 
can remain at the point of impact (partial rebound; 
Yarin, 2006). Shatter can take many different forms, 
for example prompt splash, corona splash and reced- 


ing breakup ( 

Yarin 

|2006 

Marengo et al. 

2011 

Mor- 

eira and Moita 2011 

). Additional impaction behaviours 


can involve droplets sliding or rolling along the surface 
(Mayo et al. 2015 Veremieiev et al. |2014 ), or the for¬ 
mation of rivulets, wherein the droplet slides downslope 


while spreading but recoil does not occur (Antonini et 
ah] 2014 Dong et al. 2013 2014). However, for the 


purpose of this paper we will use the three outcomes as 


outlined in Mercer et al. (2010). 


When droplets bounce or shatter on impact with 
a target, the retention of the spray at the position of 
impact is negligible compared to adhesion on impact. 
However, these bounce or shatter daughter droplets can 
redistribute to other parts of the target or nearby tar¬ 


gets (Dorr et al. 2014). Alternatively, droplets may 


bounce or shatter from non-target surfaces onto the tar¬ 
get. Thus, each impact event may positively contribute 
to retention, and it is important to understand when 
and how each takes place. 

There are several key stages common to any sin¬ 


gle droplet impaction (Attane et al. 2007 Mao et al. 


1997 Mercer et ah, 2010). First, the drop impinges on 


the surface and inertial forces are redirected from the 
impact trajectory to directions at a tangent to the sur¬ 
face. This initiates a spreading stage where the drop’s 
surface area increases, but kinetic energy is dissipated 


due to viscous forces within the fluid and at the ad¬ 
vancing contact line. If the drop has not shattered in 
response to the large inertial forces, spread continues 
until kinetic energy is depleted entirely. In some cases 
this may be the end of the droplet’s journey, with ad¬ 
hesion being the end result. In other cases, particularly 
when the surface is difficult to wet, the increased sur¬ 
face energy due to spread can cause recoil. If the drop 
recoils with enough energy, it will rebound from the 
surface (partially or completely; Rioboo et al. 2001). 
Alternatively, the drop will undergo several oscillations 
of spread and recoil but ultimately settle down to an 
equilibrium position and adhere to the surface. 

A common approach to modelling spread and recoil 
is to consider the changes in kinetic and surface energy 
as the drop oscillates on the surface. Several authors 
have developed ordinary differential equation (ODE) 
models based on energy equations, which describe how 


the radius of a drop varies during spread and recoil (At 


tane et al. 2007 Kim and Chun, 2001). Most notably, 


Attane et al. presented a semi-empirical model based on 
an assumed rimmed-disk geometry which closely pre¬ 
dicted the maximum spread factor for a range of fluid 
and impact parameters. This model was then extended 


by Mercer et al. (2010) to incorporate the effect of dif¬ 


fering contact angles for advancing and receding mo¬ 
tions, and also to include a predictor for bounce. 

Such ODE models are to be contrasted with the 
purely algebraic energy balance model by |Mao et ah] 
(1997), which not only closely predicts maximum spread 


factor (when compared to various data from the liter¬ 
ature; Mao et al. 1997), but also provides a bounce 
criterion, without the need to solve a nonlinear ODE 
for each impaction. The reduced computational time 
of an algebraic bounce model is highly advantageous 
for large scale applications where thousands of drops 
may be present at the same time, such as in multi- 
component simulations of agricultural spray scenarios 


(Dorr et al. 2014). While there are models in the litera¬ 


ture which can accurately simulate the entire impaction 


process in three dimensions (Malgarinos et ah, 2014), 


they are very computationally intensive even for a sin¬ 
gle droplet. It is for this reason that the present work 
focuses on purely algebraic energy balance models. We 
employ Mao et al.’s model to describe and predict drop¬ 
let bounce, with simple modifications to allow for the 
calculation of rebound velocity, and to apply the model 
to the case of an oblique impaction. 

The mechanisms behind droplet shatter are not yet 


completely understood (Rein and Delplanque 2008). 


As a result, much of the literature consists of empiri¬ 
cal models for the prediction of shatter based on one 
theoretical foundation in hydrodynamics: droplet shat- 
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ter occurs when the inertial forces of impact overcome 


the capillary effects of the fluid (Moreira and Moita 


20101. Empirical shatter criteria are typically simple 


equations involving dimensionless quantities such as the 
Weber, Reynolds and Ohnesorge numbers, which in¬ 


dicate whether shatter occurs or not (Mundo et al. 


1995 


Vander Wal et ah] |2006| [Rein and Delplanque 


20081. In addition to the hydrodynamic argument on 


which these criteria hinge, there is an ongoing inter¬ 
est in the literature in ideas such as the role that air 


entrapment ( 

Rein and Delplanque 

00 

o 

o 

Mandre and 

Brenner 2012 Kolinski et al. 2014) and atmospheric 

pressure (Xu et al. 2005 Riboux and Gordillo| 

2014) 


play in the onset of droplet shatter. Here we simply use 
the Mundo et al. criterion to predict shatter for normal 


and oblique impactions (Mundo et al. 19951. 


There are few studies which model the post-shatter 
behaviour deterministically (number, velocity and tra¬ 
jectory of secondary drops), although it can be achieved 
to some extent through energy balance models very sim¬ 


ilar to those used for spread and bounce (Yoon and 


Desjardin} 2006). Alternatively, empirically-based sta¬ 
tistical models may be used, such as those by |Saint-Jean| 
et al. (2004); Gigot et al. (2014) for the rain-splash dis¬ 


persal of pathogens. We do not consider post-impaction 
events in this study, as these processes are difficult to 
observe and measure accurately within the current ex¬ 
perimental framework. 

High speed video images have been used by research¬ 
ers to study the impaction process, although these have 
typically been on uniform, well defined surfaces such as 


glass 

Roux and Cooper-White 2004) 

parafilm (Wang 

et al. 2009) 

or PTFE coated microscope slide ( 

Massi- 

non and Lebeau 2012). 

Reichard et al. 

(1986), Reichard 

et al. 

(1998) 

, Fox et al. 

(1992), and Wirth et al. 

(19911 


used high-speed cameras to record droplet impaction 
on leaves. The speed and image resolution of these cam¬ 
eras, however, could not capture the dynamic impaction 
process, and the properties of daughter droplets could 
not be determined because many droplets overlapped 
within the camera target focus area. More recently. 


Dong et al. (2013) and Dong et al. (2014) developed 


a system that used two high-speed digital cameras and 
3D software to study the impact of a range of drop¬ 
let sizes and impact velocities on Hydrangea quercifolia, 
Dracaena, and Euphorbia pulcherrima, and was used to 
qualitatively describe the impaction process and mea¬ 
sure the spread of droplets on impact. 


In other recent research, Massinon and Lebeau (2012) 


and Zwertvaegher et al. (2014) studied droplets of both 


water and surfactant mixtures impacting on horizontal, 
synthetic, hydrophobic polytetrafluoroethylene. Teflon® 
(PTFE) coated slides. Drops were delivered to the sur¬ 


face with a moving agricultural nozzle, and impaction 
details (droplet size and velocity) and outcomes (ad¬ 
hesion, bounce and shatter) were extracted from high 
speed video images. Their results show a zone of tran¬ 
sition between each impaction outcome and not a clear 
boundary between them. A ‘Weber number of transi¬ 
tion’ was used to delineate the boundary between adhe¬ 
sion and bounce, and bounce and shatter. These transi¬ 
tion values were obtained from the intersection between 
the Weber number probability density distributions of 
their data. In some cases (when surfactant concentra¬ 
tion was sufficiently increased), bounce outcomes were 
not observed at all, and so only one Weber number of 
transition was required, describing a boundary between 


adhesion and shatter zones (Massinon and Lebeau 20121 


Likewise, Massinon et al. (2014) obtained similar re¬ 


sults for water and surfactant mixtures on PTFE and 
blackgrass leaves of various orientations. Rather than 
using Weber numbers of transition to delineate adhesion- 
bounce and bounce-shatter zones, they reported im¬ 
pact outcome probability as a function of eleven energy 
classes, where each class boundary corresponds to a 
constant Weber number. As with the Weber numbers of 


transition ( 

Massinon and Lebeau 

2012 

Zwertvaegher 

et al. 

2014 

1, 

Massinon et al. 

(20D 

;)’s impact outcome 


probabilities in each energy class were experimentally- 
determined. Spraying distilled water was found to result 
in less adhesion, more bounce, less shatter in Cassie- 
Baxter regime and more shatter in Wenzel regime than 
when a spray mixture of 0.1% Break-Thru S240® was 
applied. For impaction on blackgrass leaves, |Massinon| 
et al. (20141 found that the leaf roughness pattern clearly 


affects drop impact behavior compared with an artifi¬ 
cial surface. From these results it seems reasonable to 
expect that impaction outcomes may be more variable 
on natural leaf surfaces due to their inherent variability 


Nairn et al. 2014 

Forster et al. 

2005 

Taylor 2011 


Experimental techniques such as those outlined by 


Massinon and Lebeau 

2012 

, Zwertvaegher et al. 

(2014) 

Massinon et al. (2014 

and 

Dong et al. (2014) demon- 


strate a difference in impaction outcomes between dif¬ 
ferent surfaces and spray mixtures. However, it is desir¬ 
able for future modelling of whole plant spraying, and 
also pathogen spreading scenarios, to be able to predict 
these effects without being totally reliant on experimen¬ 
tal data for each possible application scenario. 


In this paper we report on experimental data taken 
from high speed video images of droplets impacting on 
cotton, wheat and rice leaves. In an attempt to more 
closely mimic conditions in the field, we use three dif¬ 
ferent commercial agrichemical spray nozzles and four 
different formulation types (one of which is simply wa¬ 
ter, while the other three include commonly used ad- 







































































































































































































4 


Gary J. Dorr et al. 


Track sprayer with single nozzle (moving into page) 



A 


Cover 


10 ^ ^ 

o 

o 

LT) 



Fig. 1 Schematic layout of high speed video experiments 


juvants). The experimental details are provided in Sec¬ 
tion [21 For both horizontal leaves and leaves held at 
45°, the outcome (adhere, bounce or shatter) for each 
droplet impaction was recorded, as well as the droplet 
size and speed just before impaction. The mathematical 
approaches outlined above (and described in Section]^ 
are tested by plotting bounce and shatter boundaries 
as a function of droplet size and speed, and compar¬ 
ing with the experimental results (Section]^. Section]^ 
includes a discussion of our work. We find broad agree¬ 
ment between the data and predictions, but also high¬ 
light cases in which observations are not backed up by 
the models as they stand. This research provides confi¬ 
dence that simple and computationally efficient math¬ 
ematical models can be used to investigate retention of 
pesticides on whole plant from typical pesticide applica¬ 
tions (Dorr et al., 2008 2014), but also points to aspects 
of the modelling which needs further refinement. 


2 High speed video analysis 

2.1 Experimental set up 

Droplets with a range of diameters and velocities were 
recorded impacting on three plant species with a high 
speed video system at a frame rate of 5000 f/s and shut¬ 
ter speed of 1/20000 s. The system consisted of a Light¬ 
ning RDT high-speed camera (DRS Technologies, 
Inc., USA), an HL-250 high-intensity discharge light 
with maximum power of 1, 250W (Shanghai Anma In¬ 
dustry Co., Ltd, China) and integrated with MiDAS 
2.0 software modules (Xcitex, Inc). The experimental 
layout is shown in Fig[^ 

During each test run, a target leaf was fixed on a 
sample board and placed on an aluminum alloy frame 
that enabled the angle of the leaf to be altered. The 


camera focused on the leaf surface to capture the im¬ 
paction process. To protect the camera and light from 
spray and to limit the number of out of focus droplets, 
a plastic board with a slot 0.01 m wide and 0.1 m long 
was used as a cover. The slot was positioned over the 
leaf to allow droplets to go through and impact on the 
target at the camera focal plane. 

By using a slot to limit the number of droplets, to¬ 
gether with a low application rate that is typical of field 
applications, we were able to ensure that the majority 
of the impacting droplets hit a dry area on the leaf. To 
verify that most of the droplets impact an area of the 
leaf where there was no previous impact, we sprayed 
water under the same operating conditions on water 
sensitive cards (75 x 25 mm) that change color from 
yellow to blue when a droplet impacts on their surface. 

A figure demonstrating these results is included in the 
supplementary material. 

After recording the spray, the movement of droplets 
that were in good focus prior to impact was manu¬ 
ally tracked using the MiDAS software and the im¬ 
pact velocity of the droplet was calculated. It was then 
noted from the recorded images if that droplet adhered, 
bounced or shattered on impact. To measure the drop¬ 
let size, a public domain image processing program (Im¬ 
aged 1.48c, developed at the National Institutes of Health) 
was used. Towards the end of the frame sequence, droplets 
that were suspected to have hit a wet area (due to a vis¬ 
ible change in behavior of the impacting droplet) were 
excluded from the analysis. 


2.2 Spray methods and formulations 


The distance from the nozzle tip to the target was 0.5 
m. The spray pressure was 0.3 MPa, supplied by a 
TYW-2 air compressor (Suzhou Tongyi Electromechan¬ 
ical Co., Ltd., China). The droplets were produced by 
three nozzle types (IDK120-02, STllO-02 and TR80- 
02 nozzles. Techier GmbH, Germany) whose respec¬ 
tive Volume Median Diameters (VMD) with water were 
273 ^m, 174 /im and 171 /rm, respectively, measured 
by a SympatecHelios laser-diffraction particle-size anal¬ 
yser (Spraytec laser-diffraction particle-size analyser, 
Malvern Instruments Ltd., UK) at the stated spray dis¬ 
tance and pressure ( Wang et al.[ 20141. The nozzle was 
mounted on a track sprayer, operated at a constant 
speed of 1.0 m/s over the test zone. 

The following solutions were tested: 


— Water 

— 0.05% NEIOO (Spray adjuvant, Numen (Beijing) In¬ 
ternational Biotech Go., Ltd., China) -|- water 

— 0.1% Tween 20 -I- water 
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— 0.05% Silwet 408 (Spray adjuvant, Numen (Beijing) 
International Biotech Co., Ltd., China) + water 

Tween 20 and NFIOO are conventional nonionic surfac¬ 
tants. Silwet 408 is an example of what is often referred 
to as a superspreader due to its very low surface ten¬ 
sion and ability to spread over the leaf surface after 


1991 

Stevens et al. 

1993 

Forster et al. 

2005 


impact (Ananthapadmanabhan et al. 

1990 

Nikolov et 

al. 

2002| Venzmer, 2011 Wang et al. 

2013 

Zhang et 

al. 

2006). 


Each solution was sprayed onto leaves from three 
crops (cotton, rice and wheat) for both horizontal leaves 
(0°) and leaves at a 45° angle, to give 72 treatments. 
Plants were grown in a glasshouse for approximately 
one month prior to tests. Every treatment was repeated 
twice. The data from the three nozzle types was com¬ 
bined so that impaction results for each leaf type, spray 
mixture, and leaf angle combination were obtained from 
six separate leaves, each from a different plant. 

The contact angle for each of the spray mixes on 
cotton leaves was initially measured by placing a 1 /xL 
droplet with a micro syringe on a section of a cotton 
leaf, viewing side on through a microscope, manually 
measuring the height and base width of the droplet, and 
calculating the contact angle by the method in |Mack| 
(1936). It was found that by the time the leaf section 


was positioned on the microscope, the microscope was 
focused and the measurements taken, droplets on wheat 
and rice leaves had flattened, resulting in lower than 
anticipated contact angle measurements. To overcome 
this issue, the high speed video camera was then used to 
record images of the droplet as it was applied to the leaf. 
The contact angle was determined from the resulting 
images by using the low-bond axisymmetric drop shape 


analysis (LBADSA) plug-in to Imaged (Stalder et al. 


2010 ). 


To delineate impaction thresholds, a value for ATcrit, 


as described in § 3.1 was calculated using contact an¬ 
gles of either 20% or 50% aqueous acetone solution via 
the method in Forster et al. ( 2010[ ). The ATcrit value for 
each plant type is shown in Table 

Dynamic surface tension (DST) was measured with 
a Bubble pressure Tensiometer-BP2 (Kriiss, Germany). 
It is well established that DST, rather than equilibrium 


surface tension, is better correlated with retention (An¬ 


derson and Hall 19891. However, the most appropriate 


time scale for the DST measurement is not yet fully 
resolved. On one hand it may be argued that the drop¬ 
let surface is instantaneously ‘fresh’ at impact due to a 


significant deformation of the fluid free surface (Ander¬ 


son and Hall 1989), and so the most appropriate time 
scale may as early as 10 ms. Others, however, argue 
that the DST measurement should correspond to the 


the most appropriate time scale for DST has not been 
resolved, we consider values at a surface age of both 10 
ms and 50 ms. The former represents a fresh droplet 
surface, while the latter is the average time between 
atomization and impact. Since our DST data began at 
roughly 30 ms, the 10 ms estimate was obtained by ex¬ 


trapolation from a fitting of Hua and Rosen (1988)’s 


Equation (2) for the DST of a surfactant mixture. To 
arrive at the 50 ms estimate, we assume the average 
velocity of a droplet from nozzle to target is roughly 
10 m/s (all of the measured impact velocities in Figs[^ 
and are less than 10 m/s, but drop velocities close to 
the nozzle are greater than 10 m/s, Wang et al., 2015), 
and the distance from nozzle to target is 0.5 m. 

For comparison, we have also considered DST at 
900 ms, which is essentially the equilibrium value for 
our three spray mixtures. The validity of using DST at 
these three ages is discussed in Sections |5.1| and |5.2| 
The surface tension values for each spray mixture are 
shown in Table [T] 


3 Droplet impaction models 


In Dorr et al. (2014), simple mathematical models were 


presented to predict the outcome of a droplet impact¬ 
ing a horizontal leaf surface from directly above (nor¬ 
mal impaction). The three possible outcomes consid¬ 
ered were shatter, bounce, and adhesion. The shatter 
criterion (based on Mundo et al. 1995) is summarised 
below in § |3.1[ and in § |3.2[ an extension to the bounce 
criterion (based on Mao et al. 1997) is provided to allow 
for different angles of the leaf surface and the droplet 
impact trajectory (oblique impaction). The key param¬ 
eters describing droplet impaction are the impact ve¬ 
locity (E), droplet diameter before impact (D), fluid 
density (p), dynamic fluid viscosity (p), fluid-air sur¬ 
face tension (ct), and equilibrium contact angle (0e)- 
Note that each spray mixture is assumed to have the 
density and viscosity of water (p = 1000 kg/m^ and 
/i = 9.78 X 10“^ Pa-s); it is only surface tension that 
changes with the addition of an adjuvant (see Table 
for specific surface tension values). 


3.1 Shatter criterion 

The process of droplet shatter can be quite variable, 
depending on the nature of the fluid and the surface 
being impacted. Types of shatter can be categorised as 


prompt, crown, and receding, to name a few (Yarin 


time to impact (Anderson and Hall, 1989 Wirth et al. 2006 Marengo et al. 2011). This wide variety, together 
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Table 1 Surface tension and contact angles of the test spray mixtures on leaf surfaces. Standard deviation from 3 replicate 
measurements shown in brackets. 



Surface tension a (N/i 

m) 


Contact angle 

Mix 

10 ms 

50 ms 

900 ms 

Cotton 

Wheat 

Rice 

Water 

0.05% NFIOO 
0.10% Tween 20 
0.05% Silwet 408 

0.0716 (0.0001) 
0.0640 (0.0002) 
0.0457 (0.0001) 
0.0514 (0.0002) 

0.0716 (0.0001) 
0.0569 (0.0002) 
0.0374 (0.0001) 
0.0359 (0.0002) 

0.0716 (0.0001) 
0.0482 (0.0002) 
0.0320 (0.0001) 
0.0194 (0.0002) 

54° (10) 
48° (3) 

55° (7) 

* 

132° (2) 
110° (4) 
116° (4) 
69° (9) 

129° (6) 
126° (3) 
108° (4) 
73° (2) 



A'crit (dimensionless) 

150 

69 

74 


*Complete spread so value of 1° used. 


with the fact that the physical processes occurring dur¬ 
ing shatter are not entirely understood, makes it a very 
difficult subject to accurately model and predict using 
simple models. 

It is widely agreed that a prompt splash occurs when 
the inertial forces created by impact are sufficient to 
overcome the capillary forces which hold the fluid to¬ 


gether (Moreira and Moita 2010). This takes place dur¬ 


ing the spreading process, immediately after impaction. 
Several authors have developed criteria to predict shat¬ 


ter based on this theory, such as Mundo et al. (1995) 
who proposed the relation 


K = , 


( 1 ) 


where Wen = R,en = pV^D/^ are the 

dimensionless Weber and Reynolds numbers computed 
with the component of velocity normal to the impacted 
surface 


V^ = V sin a. 


( 2 ) 


Here a is the angle between the leaf surface and the 
incoming trajectory of impact (0 < a < 90°). Thus 0 
is valid for both normal and oblique impactions. 

A droplet is predicted to shatter on impact if 


K>K, 


( 3 ) 


where ATcrit is a critical value related to the properties 
of the surface being impacted. Mundo et al. found that 
A'crit = 57.7 correlated well with their data (which in¬ 
volved impactions of three different fluids onto smooth 
and rough surfaces, which rotated to mimic impact at 
an angle); however, this value could not be expected to 
translate to the data presented in this study. |Mao et| 
al. ( 1997|, for example, demonstrated that ATcrit could 


be as large as 152 for particular surfaces. Therefore cri¬ 
teria like 0 typically require that laborious adhesion 
and shatter experiments are performed in order to fit 
a suitable value of A'crit to any new data set. To re¬ 
move this limitation, 


simple method for estimating ATcrit based on contact 
angle measurements of standardised formulations. In 
this way, the value of ATcrit adapts to reflect the wet¬ 
tability properties (including roughness effects) of the 
surface, while AT describes the fluid and droplet prop¬ 
erties. Hence the novelty in this study is in the use of 
Forster et al. ’s method to extend upon [Mundo et al.| 
(|1995|)’s criterion. 


Note that as the angle of impact a decreases, 14 will 
in turn decrease, leading to a smaller calculated value of 
K. The implication of this trend is that shatter becomes 
less likely with a smaller impact angle. 

It is worth noting that ([^ may be equivalently writ¬ 
ten as K = OhRe^'^^, where Oh = fij{pgD)^^'^ is the 
Ohnesorge number. Further, 0 may be written as Oh > 
ATcritThus, as summarised by Rein and Delplanque 


(20081, this criterion belongs to a family of models pre¬ 
dicting that shatter occurs when Oh > a Ref, with a 
constant and f3 = —BjA. Other models in this family 


are due to Cossali et al. (1997) and Yarin and Weiss 


(1995). In these works, the focus was on droplet im¬ 


paction on a prewetted surface. Our experimental set 
up is such that droplets usually impact dry areas, al¬ 
though towards the end of the run there is an increased 
chance of impacting on portions of the leaf surface that 
had previously been wet by other droplets. Thus we do 


and Weiss (1995). 


not adopt the models of Cossali et al. (1997) or Yarin 


Mundo et al.’s shatter criterion, provides a method 
for predicting whether splash will occur or not. How¬ 
ever, it does not begin to address the secondary distri¬ 
bution of the spray, such as the number of secondary 
droplets produced in the splash, and their individual ve¬ 
locities and trajectories from the impaction site. Such 
post-shatter information is not considered in the cur¬ 
rent study, but relevant models have been published in 


Forster et al. 

(2010 

) devised a 

Yoon and Desjardin 

(2006 

) and 

Dorr et al. 

(2014 
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3.2 Bounce criterion 


In the event that (§ predicts that shatter does not 
occur, a secondary criterion is required to differentiate 
between the possible subsequent outcomes of bounce 
and adhesion. Our bounce criterion hinges on work by 


Mao et al. (1997), with a few modifications. 


Mao et al.’s bounce model considers several key stages 
of the impaction process: 


(a) Before impact. A spherical droplet of diameter D 
is travelling at velocity V towards a solid surface. 

(b) Maximum spread. After impact, the droplet spreads 
radially to a maximum diameter. A cylindrical disk 
geometry is assumed. 

(c) Maximum recoil. After maximum spread, the drop¬ 
let recoils toward the point of impact. 

(d) Equilibrium. If the droplet does not bounce, it will 
settle to an equilibrium formation on the surface. 


In addition to these four, a fictitious stage is introduced 
for the purpose of predicting bounce: 

(e) Bounce criterion. A spherical droplet momentarily 
at rest above the surface. 


Equating stages (a) and (b) leads to a cubic polyno¬ 
mial describing the maximum spread diameter of the 
droplet. It is assumed that, if bounce occurs, it takes 
place after stage (c), and instead of (d), in response to 
a sufficiently rapid recoil. A comparison between stage 
(c) and the fictitious stage (e) determines whether de¬ 
tachment from the surface (bounce) is the energetically 
favourable outcome after recoil. A successful bounce is 
indicated by a positive value of an ‘excess rebound en¬ 
ergy’: 


Aere > 0. 


(4) 


If Aere < 0, then the droplet is predicted to adhere to 
the substrate. 


Dorr et al. (2014| expanded upon Mao et al.’s model 


by using Eere to calculate the exit velocity of the 
bouncing droplet. However, the model fundamentally 
describes the spread and bounce of a drop impinging a 
solid horizontal substrate from directly above (normal 
impaction). This is a major limitation in a spray sce¬ 
nario where leaf surfaces may have any orientation, and 
are impacted by spray droplets from various directions. 
Here we present a modification of the bounce criterion 
to account for oblique impaction. 

The impact angle parameter a, is used to quantify 
the ‘obliqueness’ of an impact, with a = 90° for a nor¬ 
mal impaction and smaller a being more oblique. The 
orientation of the leaf surface itself is not taken into 
account explicitly because it is assumed that the mag¬ 
nitude of the inertial forces acting on the droplet during 


impaction, over the short timescale of spread and recoil, 
is much greater than the effect of gravity on the droplet. 
This is reasonable when considering typical spray-sized 
droplets of very small mass. 

The modelling for oblique impact is kept as simi¬ 
lar to the case of normal impaction as possible. Indeed, 
the criterion outlined here reduces to that in IDorr et al.l 


(20141 when a = 90°. While a cylindrical disk geometry 


is assumed by Mao et al. (1997) for a normal impaction. 


an elliptical disk geometry is now taken for the spread¬ 
ing drop, with the major axis in the impact direction. 
There is some probability of an oblique impaction lead¬ 
ing to rivulet-type stretching, sliding, or rolling down¬ 
hill (particularly for small a and highly hydrophobic 
surfaces; Antonini et al. 2014). However, it has been ob¬ 
served that a drop impacting an oblique surface can also 
exhibit the conventional spread and recoil behaviour 


that we assume here (Liu et al. 20141. Further, ellipses 


have been used in the literature to approximate the 
shape of oblique spreading patterns in a few applica- 


tions ( 

Kang and Ng 

2006 

Chen and Wang 

2005 

Zen 

et al. 

2010). Thus, we acknowledge that the use of an 


elliptical geometry and the assumptions behind it is an 
ad hoc approach which meets our requirements of a sim¬ 
plistic criterion with realistic bounce trends. There may 
be other simple alternatives that work better, especially 
in situations where a is expected to be small. 

To quantify the amount of spread along the minor 
and major axes of the assumed ellipse, a few simple 
assumptions are made about the spreading process: 


1. The sideways spread of the droplet (in the direction 
of the minor axis of the ellipse, perpendicular to the 
impaction direction) can be quantified using the nor¬ 
mal component of impact velocity. 


Although the majority of spread will be in the impact 
direction, some spread will occur in the sideways di¬ 
rection as well. The sideways spread will decrease as a 
decreases, along with the normal component of impact 
velocity, 14. 


2. The area of the droplet at maximum elliptical spread 
is equal to what it would be if impacting at a = 90° 
with the same impact velocity. 


The principle here is that the droplet should undergo 
roughly the same extent of deformation when impact¬ 
ing a surface at velocity H, regardless of whether the 
impaction is oblique or normal. 


3. The retraction of the droplet after maximum spread 
will depend only on the major axis. 

A droplet bounces if there is sufficient surface energy 
to cause the droplet to recoil rapidly. If the spread is 
elliptical, the major axis will determine the likelihood 
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of rapid recoil while the minor ellipse will have little 
contribution to the process. 

The method employed for predicting oblique bounce 
is as follows. First, the impaction is treated as if it is 


normal {a = 90°) (Mao et al. 1997 Dorr et al. 2014). 
The maximum spread diameter Dnormai is found via the 
cubic polynomial 

f Tlnormal \ 


1 We°®^ 

-(1 — COS^e) + 0.2 J^^o.33 


\ D J 


/We ^ 

l2 


/ ^normal A 

J 




( 5 ) 


which can be solved exactly for Dnormai- This corre¬ 
sponds to Equation (17) of (Mao et al. 1997), and 


makes use of a stagnation point flow analogy to empir¬ 
ically quantify viscous dissipation during the spreading 
process. 

If a real solution to ([^ does not exist, or Dnormai 
is calculated to be less than D, we set Dnormai = D as 
a minimum bound. Once D„ormai is known, the second 
assumption implies that the spread area of the droplet 
in the oblique impaction is 


Area = —D 


normal * 


( 6 ) 


Using the first assumption, the extent of sideways spread 
of the elliptical droplet is calculated with the normal 
component of the velocity (Equation leading to 
another cubic equation for maximum spread: 

3 


1 We°®^'' 

-(1 — COS0e) + 0.2 j^^o.33 


Dn 


Wen 

12 


+ 1 


D„ 


D 


D 




( 7 ) 


to be solved for Dminor under the condition. The no¬ 
tation Wen and Ren denotes the Weber and Reynolds 
numbers computed with Vn rather than V. If Dminor < 
D, then we set Dminor = D as a correction. 

Now, the area of an ellipse is given by 


Area — ^ DmlnorTlmajor ; 


( 8 ) 


and so the second assumption coupled with §> allows 
the calculation of D„ 


as 


D, 


major 


n2 

normal 


( 9 ) 


Using the third assumption, the excess rebound energy 
of the droplet is calculated by 


Eere — 


4 major 


2 

{1 - cosOe) + -TT -— 


major 


— 0.127rD^(T 

— iraD^ . 




D 


2.3 


(I — COSde) 


0.63 


( 10 ) 


Bounce occurs if the criterion ([^ is satisfied (otherwise 
adhesion is predicted). Equations ([^-(|l0|) incorporate 
properties of the fluid (density, viscosity, and surface 
tension) and droplet (velocity and diameter) through 
the Weber and Reynolds numbers. Meanwhile, the equi¬ 
librium contact angle 9e describes the energetics of the 


liquid-solid-gas system. Mao et al. (1997) found that 


9e was sufficient to account for the effects of surface 
roughness without explicitly including a roughness pa¬ 
rameter. The exit velocity of the droplet is 


14xit — 


I12E, 


ERE 


npD^ 


( 11 ) 


The direction that the droplet bounces is assumed to 
simply be a mirror of its incoming trajectory. The va¬ 
lidity of this assumption is discussed in § |^ 

The impaction models described here are best vi¬ 
sualised by a bounce or shatter boundary. Eor a given 
plant and spray formulation, one could run the bounce 
criterion for an appropriate range of V and D, and 
then on the V versus D plane delineate regions of pre¬ 
dicted adhesion and bounce with a curve representing 
the boundary between the two outcomes. This bound¬ 
ary corresponds to Eere = 0. Similarly, K = ATcrit 
leads to a shatter boundary which may be plotted on 
the same plane. The interpretation of the bounce (or 
shatter) boundary is that a spray droplet with initial 
velocity and diameter falling above the curve results in 
a positive prediction for bounce (or shatter). Recall that 
bounce may only occur if there was no shatter during 
the droplet’s spreading process, and so shatter always 
takes precedence over bounce in the impaction model. 


4 Results 

Experimental impact outcomes for droplets on cotton 
and wheat leaves are shown in Eigj^for horizontal leaf 
surfaces, and Figfor leaf surfaces inclined at an angle 
of 45°. For each leaf type, the results for the four spray 
mixtures (water. Tween 20, Silwet 408 and NF 100) are 
shown separately. Each of the subplots combines data 
from the three nozzle types, and on average 75 impact 
outcomes, on six individual leaves, from six different 
plants, are recorded in each. Additional data for rice is 
provided in the online supplementary material. 

Each data point in Figs and [^represents a single 
droplet impaction, plotted according to the measured 
droplet diameter and velocity on impact. Three symbols 
are used to indicate the observed impaction outcomes: 
a (red) square for droplet shatter on impact, a (blue) 
circle for droplet bounce on impact, and a (black) cross 
for droplet adhesion on impact. All three outcomes are 
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Cotton - horizontal 


Wheat - horizontal 


water 


water 





Tween 20 


Tween 20 



Silwet 408 


Silwet 408 


a 

Q 



V (m/s) 


V (m/s) 


■ Shatter • Bounce x Adhere - Bounce boundary - Shatter boundary 

Fig. 2 Impaction observations for horizontal wheat and cotton leaves and a variety of formulations. The corresponding 
predictions are illustrated with bounce (i?ERE = 0) and shatter {K = K^rit) curves. The top shatter curve in each subplot 
corresponds to DST at a surface age of 10 ms, the middle curve to 50 ms, and the bottom curve to 900 ms. The bounce curves 
are computed with DST at 50 ms. 
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Cotton - 45° angle 


Wheat - 45° angle 


water 


water 





Silwet 408 


Silwet 408 



■ Shatter • Bounce X Adhere - Bounce boundary - Shatter boundary 

Fig. 3 Impaction observations for wheat and cotton leaves inclined at 45°, and a variety of formulations. The corresponding 
predictions are illustrated with bounce (i?ERE = 0) and shatter {K = K^rit) curves. The top shatter curve in each subplot 
corresponds to DST at a surface age of 10 ms, the middle curve to 50 ms, and the bottom curve to 900 ms. The bounce curves 
are computed with DST at 50 ms. 
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well represented in the results for wheat leaves, while 
the results for cotton leaves show a majority of adhesion 
observations with relatively few shatter and bounce ob¬ 
servations. For all leaf and spray mixture types, there 
is not always a clear grouping of the three types of im¬ 
paction outcome. This is particularly true for bounce 
and adhesion observations, where the data are often 
mixed in the lower left region of the plots (correspond¬ 
ing to smaller diameters and velocities). The shatter 
observations, on the other hand, are generally grouped 
in the upper right region (corresponding to larger di¬ 
ameters and velocities) with little overlap with bounce 
and adhesion observations. It is interesting to note that 
there is not a clear difference in the distribution of shat¬ 
ter, bounce and adhesion results between the horizontal 
and 45° leaves. 


For comparison of the experimental data with im¬ 
paction criteria, shatter boundaries {K = Kent, from 
the model in § 3.11 are indicated by a (red) dashed 
line, and bounce boundaries (I^ere = 0, as described 
in § 3.2) by a (blue) solid line. Recall that shatter is pre¬ 
dicted when a droplet’s impact diameter and velocity 
fall above the shatter boundary, so that is satis¬ 
fied. On the other hand, bounce is predicted when the 
droplet’s attributes place it between the bounce and 
shatter boundaries, because shatter takes precedence 
over bounce. In those plots within Figs and where 
only a shatter boundary is depicted, the bounce bound¬ 
ary has been omitted because it fell above the shatter 
boundary. Droplet adhesion is predicted when the di¬ 
ameter and velocity information falls below or to the 
left of a bounce boundary, or below and to the left of the 
shatter boundary for cases without a bounce boundary. 
Three shatter boundaries are shown in each subplot of 
Figs and these correspond to DST measurements 
at 10 ms (top shatter boundary), 50 ms (middle shat¬ 
ter boundary), and 900 ms (bottom shatter boundary). 
Bounce boundaries are shown for a DST measurement 
at 50 ms only. This is for two reasons. First, a change in 
DST made little visible difference to the position of the 
bounce boundary, and so does not affect interpretation 
of the results. Second, the most appropriate DST mea¬ 
surement is even more uncertain for the bounce process 
than it is for shatter, due to the larger timescale over 
which it occurs and the potential for DST to vary sig¬ 
nificantly over that time. 


Note that there are a few important trends in the 
positions of the shatter and bounce boundaries in Figs[^ 
and[^ First, the position of the shatter boundary low¬ 
ers as the surface tension of the spray mix is decreased. 
This is true both between subplots (when comparing 
spray formulations), and within subplots (when com¬ 
paring DST measurements). This means that the prob¬ 


ability of shatter increases with decreasing surface ten¬ 
sion, because the droplets have a lower surface energy 
and are therefore more prone to break-up. Second, the 
position of the shatter boundary rises as the leaf sur¬ 
face is tilted from the horizontal to 45°. This means 
that shatter becomes less probable in an oblique im¬ 
paction. Finally, the bounce boundary slightly lowers 
as the leaf is tilted to 45°, indicating a higher probabil¬ 
ity of bounce in an oblique impaction. 

In Fig|^a), the observed angle of bounce, as mea¬ 
sured from the nominal surface normal, is contrasted 
with the angle of impact. The (red) line running diag¬ 
onally through the plot indicates the expected bounce 
angle if rebound was to occur as a mirror image of the 
incoming trajectory. The distribution of the majority 
of data points below the line indicates that the bounce 
angle tends to be greater than the impact angle, im¬ 
plying that the assumption made in § |3.2| about these 
angles being equal is not strictly applicable. This trend 
can be attributed to energy loss in the impact process 
(spread and retraction of the droplets). It should be 
noted that while these results are relative to the nom¬ 
inal surface normal (normal to either a horizontal or 
45° surface), the leaf surfaces have irregularities that 
would cause the actual surface normal at the point of 
impact to vary from the nominal surface normal. This 
difference would contribute to the scatter in the data, 
and may account for the cases where the bounce angle 
is less than the impact angle. 

Figurej^b) compares the measured velocity of bounc¬ 
ing droplets to the predicted velocity calculated by Equa¬ 
tion 0- The (red) line illustrates where data points 
would fall if there was perfect agreement between the 
experiments and prediction. There is a slight tendency 
for droplets with higher observed bounce velocities (> 1 
m/s) to fall below the line, indicating that bounce ve¬ 
locity may be underestimated by the model. The overall 
agreement is quite good, however, with the scatter of 
data points generally centered about the line. 


5 Discussion 

While there exist a variety of sophisticated mathemat¬ 
ical models that are designed to simulate the entire 
droplet impaction process, the majority of these are 
computationally intensive even for a single droplet. On 
the other hand, we are motivated by large scale appli¬ 
cations where thousands of drops may be present at 
the same time, such as in multicomponent simulations 
of agricultural spray scenarios. Thus we have a need 
for simple, easy-to-compute models and, as such, have 
focused on algebraic energy balance criteria and an ad 
hoc type approach. While there is scatter present in the 
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-0 
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(b) 



Impact angle (°) Measured bounce velocity (m/s) 

Fig. 4 (a) Measured impact angle and bounce angle for all droplets that were observed to bounce. The (red) line shows the 
expected angle if bounce were to occur as a mirror image to the incoming trajectory, about the surface normal, (b) Measured 
velocity of droplet after bouncing off a surface compared to the predicted velocity. The (red) line indicates where measured 
velocity would the same as the predicted velocity. 


experimental data presented here, and distinct group¬ 
ing of different impact outcomes (adhesion, shatter or 
bounce) is not always observed in Figs and we find 
the overall predictions of the simple energy balance cri¬ 
teria compare well to the measurements. 

Most of the experimental studies on droplet im¬ 
paction to date have used uniform synthetic (dry) tar¬ 
gets and mono-sized droplets. However, in this work, 
we have used real leaves, with all their inherent nat¬ 
ural variability (for example, hairs, trichomes, stom¬ 
ata, veins, wax structures), and commercial agricultural 
spray nozzles that result in a range of droplet veloci¬ 


ties, sizes and trajectories being produced (Dorr et al. 


20131. The nature of the spray experiments introduces 


the further complication that, while the leaf surfaces 
are initially dry, they do become partially wet through¬ 
out the run as drops accumulate on the surface. There 
are some measurement errors associated with the cam¬ 
era field of view and using a single camera to record a 
3D process (Dong et al. 2013 2014), although we have 
tried to minimise these by only selecting measurement 
of drops that remain in good focus throughout the im¬ 
paction process. 


5.1 Shatter 

The shatter boundary tends to work reasonably well 
for the situations tested, although there are some dis¬ 
crepancies where most shatter is observed below the 
predicted shatter line (see water on 45° cotton and 
water and NF 100 on 45° wheat in Fig [^. There is 
also little difference between shatter boundaries for the 
three DST values, apart from in the case of the su¬ 
perspreader mix where the boundaries are more spread 
apart. Therefore it is difficult to make any conclusions 


about the most appropriate timescale for DST with re¬ 
spect to shatter. The occurrence of droplet shatter in¬ 
creases with diameter and velocity, as evidenced by the 
concentration of shatter observations in the upper-right 
of the plots. This trend is more obvious on the difficult 
to wet wheat leaves than the easy to wet cotton leaves 
where shatter observations were more scarce. 


In theory we might expect to find that the number 
of shatter observations increases as the surface tension 
of the spray mix is decreased, and that a higher impact 
velocity would be required for shatter as the surface 
tension is increased. This is due to the hydrodynamic 
principle that the inertial forces of impact must over¬ 
come the capillary effects, which are weaker as the sur¬ 
face tension decreases. There is perhaps some evidence 
in Figsand that there are fewer shatter observations 
for water on wheat than the other spray mixes, and that 
those observations generally occur at higher velocities. 
However, more data is required for confirmation. 


As stated above, there is a tendency for some drops 
to shatter when their velocity and diameter fall below 
the shatter boundary. Interestingly, these discrepancies 
tend to be present for the oblique impactions (Fig 
more than the normal impactions (Fig[^. It is likely 
that at least some of these droplets hit a previously 
wet area on the leaf and are hence more prone to shat¬ 
ter (Vander Wal et al., 2006), which is supported by 
the fact that these cases often appear later in the run 
where there is a higher chance of impacting a previously 
wet area. As the camera view is looking side on to the 
leaf surface, it is often not possible to determine if the 
impact area has been wet by previous droplet impacts 
or not. 
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5.2 Bounce 

Due to the overlap of bounce and adhesion results in 
Figs and there is not always a distinct region of 
bounce only observations in the data as there was for 
shatter, but the expected trends are still present overall. 
For example, the bounce criterion predicts that droplets 
will never bounce on the more easily wettable cotton, 
or with superspreader spray mixtures such as Silwet. In 
the data, only a few droplets were observed to bounce 
for these cases, particularly when compared to the num¬ 
ber of observations for bouncing of water droplets. 

Many bounce observations can be seen below the 
bounce boundaries in Figs and implying that the 
probability of bounce is higher than was predicted by 
the bounce criterion. This is particularly clear in the 
plot for water on a horizontal wheat leaf in Fig where 
many bounce outcomes were observed for small drops 
of low impact velocity, but not predicted by the bounce 
criterion due to low kinetic energy. 

This discrepancy between the predictions and data 
could be due to the fundamental assumptions of the 


face on impact. In reality, there is an extremely small 
air layer between the drop and substrate which must 
first be drained before contact between the two can oc- 


Liu et al. 

2015 

de Ruiter et a.| 

2015 


the discrepancy between the predictions and data could 
be due, at least in part, to complex air layer dynamics 
that are not taken into account in the bounce equations 


presented in § 3.2 


Kolinski et al. (2014a) demonstrated that, on even 


a perfectly hydrophilic substrate, the air layer can pre¬ 
vent contact and cause a drop to rebound. However, this 
was only true for an atomically smooth substrate, and 
it was shown that a drop would inevitably make con¬ 
tact with a surface with defects larger than the thick¬ 
ness of the film (as would be true for most surfaces, 
particularly natural surfaces such as leaves). Air layer 
effects have also been show to enhance the rebound of 


drops impacting on liquid. For example, Couder et al. 


(20051 demonstrated bouncing on the surface of a ver¬ 


tically vibrated liquid bath. Of more relevance to our 


study, Gilet and Bush (2012) described the presence of 


bounce criterion, adapted from Mao et al. (1997). Mao 


et al. acknowledge that the precision of the model de¬ 
creases significantly for low impact velocities {V < 1 
m/s) due to the assumption of a cylindrical disk spread¬ 
ing geometry, and a stagnation point flow analogy for 
viscous dissipation. 

One consideration may arise from the bounce ob¬ 
servations common to the context of superhydropho- 


an air layer aiding the rebound of small, slow-moving 
drops, which were impinging a surface prewetted with 
a viscous fluid. They observed that a compressible air 
layer between the drop and underlying film could act as 
a cushion, promoting bounce as long as the layer was 
not penetrated. On the other hand, faster drops could 
break through the air layer, make contact with the wet¬ 
ted surface, and adhere. These observations could rep¬ 
resent impactions occurring towards the end of our ex- 


bic surfaces. In such studies (Richard and Quere 2000 

periments, when there is a small possibility of droplets 
hitting a previously wet area. 

Okumura et al.| |20031 |Liu et ai.| 2014| Moevius et al. 

2014), the substrate is typically textured to induce a 

Gilet and Bush 

2012) reported that enhanced bounce 


Cassie-Baxter wetting regime, minimising both contact 
angle hysteresis and contact area between fluid and sub¬ 
strate on impact. As a result, drops impacting such 
surfaces have a greatly enhanced tendency to bounce. 


Richard and Quere (2000) observed that bouncing drops 


with We <C 1 (kinetic energy signihcantly smaller than 
surface energy) had very little shape deformation on 
impact, and a high restitution coefficient (e = 14 xit/f^) 
due to efficient transfer of kinetic energy to surface en¬ 
ergy. Viscous dissipation within the fluid was shown 
to be negligible. The largest contact angle considered 
in this study is 132°, so our surfaces are not consid¬ 
ered superhydrophobic, but these results illustrate why 
drops with low kinetic energy may not be adequately 
described by the current model. A different geometry 
for the impacting drop (such as an oblate spheroid, 
rather than a cylinder) may be required for this regime, 
as well as a reevaluation of the dissipation term. 

Another significant assumption of the bounce crite¬ 
rion is that a drop will make contact with the leaf sur¬ 


due to the air layer effect occurred for We < 2. They 
found that when We > 15, the air layer breaks and a 
total merger occurs between the droplet and the un¬ 
derlying film. Between We = 2 and 15, a complex series 
of intermediate regimes occur. For context, curves in¬ 
dicating We = 2 and We = 15 have been added to the 
impact outcome plot for water on a horizontal wheat 
leaf (refer to Fig[^ in Fig[^ In Gilet and Bush’s obser¬ 
vations, droplets falling above the We = 15 curve were 
not susceptible to the air layer effect, and those below 
the We = 2 curve were. 

We have presented some of many possible explana¬ 
tions for the discrepancy between observed and pre¬ 
dicted frequencies of bounce. Other potential causes 
could be the presence of dust particles on the leaves 
or irregularities in the leaf surface structure, neither of 
which are accounted for in the bounce criterion. 

In terms of total spray retention on whole plants, 
this extra incidence of bounce may not have a signifi¬ 
cant influence, since these small droplets only account 
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• Bounce X Adhere - Bounce boundary 

- We = 15 . We = 2 

Fig. 5 Magnified chart of water droplets impacting a hori¬ 
zontal wheat leaf. The predicted bounce boundary presented 
in this paper is indicated, along with curves corresponding to 
We = 2,15. 


for a very small proportion of the total volume of the 
spray, and they are typically moving at low velocities 
so they are often recaptured by the same leaf anyway. 

As with shatter discussed above, the value of dy¬ 
namic surface tension we use in our model (measured 
at 50 ms) may not be the most appropriate to describe 
the bouncing process. Since bounce occurs over a con¬ 
siderably longer timescale than shatter, the concentra¬ 
tion of surfactant molecules at the fluid surface can vary 
in a complex manner during the processes of spreading 
and recoiling, which is not yet completely understood. 
Therefore, different values of a (and consequently, 9^.) 
may need to be used for each term of Equation (10) to 
reflect the true dynamic nature of the surface tension. 


5.3 Surface wettability, roughness and texture 



Fig. 6 Scanning electron microscope (SEM) photographs at 
200-fold magnification of leaves for three crops (a) wheat, (b) 
cotton, (c) rice. 


The properties contributing to the wettability of a leaf 
can be divided into two categories: physical and chemi¬ 
cal interactions between the leaf surface and the droplet 
solution (Holloway, 19701. Most studies do not attempt 
to isolate and quantify each of these properties, but 
rather indirectly infer relative values by quantitatively 
measuring behaviour that is allied to wettability, such 
as contact angles (Nairn et al., 2011[ ) (as we have done 
in this paper). 

A leaf may be difhcult-to-wet (high contact angle of 
water, 20% or 50% acetone) even though it may have a 
chemically neutral or even relatively polar surface, due 
to the surface being extremely rough (as opposed to 
being hydrophobic) (Nairn et al., 2011 2015). 

In order to demonstrate the variability of the leaf 
surfaces from each of the three crops at a small scale, 
scanning electron microscope (SEM) photographs at 
200-fold magnihcation (Yang, 2012) are presented in 


Eig|^ These images show: the surface of a wheat leaf is 
smooth and covered by a wax layer with hooked shape 
(unciform) and needle shaped (acicular) hair on the sur¬ 
face, but no papilla (Eig|^); the surface of a cotton leaf 
has no hair nor papilla (Fig[^); and the surface of rice 
leaf is rough with unciform hairs and also papilla (Eig 
[^). 


5.4 Leaf elasticity 

The impaction models presented in this study assume a 
stationary, rigid leaf surface. This is in line with the ex¬ 
perimental observations gained from high-speed video 
imaging, in which the spray droplets were very small rel¬ 
ative to the leaf surface, the leaf was fixed on a board, 
and there was no evidence of deformation of the leaf 
surfaces on impact. It is possible, however, that the 
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impaction of large drops, the bombardment of many 
small spray drops, or the presence of wind could cause 
significant movement of leaves in a real spray scenario. 
Further, particular leaves may have compliant surfaces, 
especially when impacted by a drop with significant ki¬ 
netic energy (size or velocity). 

To incorporate leaf movement in the current bounce 
and shatter criteria, a straightforward approach may be 
to replace impact velocity with the net impact veloc¬ 
ity (sum of the droplet and leaf velocities). However, 


Lee and Kim (2004) demonstrate for the case of drop¬ 


let rebound that it is the motion of the target surface 
and surrounding air which influences the impact out¬ 
come, rather than the net change in velocity. Therefore 
air motion may need to be considered, at least in the 
spreading stage where it is most influential. 

Elasticity of the leaf surface may be accounted for 
with a small reduction in droplet kinetic energy, repre¬ 
senting the kinetic energy absorbed by the surface dur¬ 
ing impact. Pepper et al. (2008) showed that substrate 
elasticity can suppress shatter for this reason, but that 
these effects cannot be described by a traditional en¬ 
ergy balance approach due to the timescale over which 


shatter occurs. Similarly, Mangili et al. (2012) and Al- 


izadeh et al. (2013) considered drop impaction on soft 


polydimethysiloxane (PDMS) substrates, and demon¬ 
strated an absorbance of kinetic energy due to substrate 
deformation, which resulted in a significant reduction of 
droplet recoil. This suggests a lower available energy for 
rebound on compliant surfaces. In agreement with these 
studies, Cho et al. (2013) saw that their fibrous elastic 
membrane could suppress both shatter and bounce. 
Complicating these analyses are the findings of|Chen| 


and Li (2010) and Chen et al. (2011), which show that 


an air film may form beneath a drop impacting a soft 
PDMS substrate. As long as the drop velocity is within 
a critical lower and upper limit, air entrapment oc¬ 
curs between the droplet and surface due to the shear¬ 
thinning property of elastomeric PDMS. As a result, 
droplet bounce was observed on soft PDMS even when 
it did not occur for any impact velocity on rigid PDMS. 
These results, along with the discussion of air films in 
the previous section, illustrate a need for a better un¬ 
derstanding of air layer effects on droplet impaction. 


5.5 Dispersal of pathogens by rain 

While the data presented in this paper were gener¬ 
ated for typical pesticide applications, the principles 
are relevant to the impaction of rain droplets onto leaf 
surfaces and the subsequent dispersal of pathogens be¬ 
tween plants. The main differences between these two 
applications are that: 1) rain droplets are purely water. 


whereas typical spray mixtures for pesticide application 
contain additions that alter characteristics such as sur¬ 
face tension; and 2) for pesticide applications, the leaf 
surfaces are mostly dry and the volumes used are usu¬ 
ally such that the surface does not become fully wet, 
while impaction during rain events are often onto fully 
wet surfaces. It would be expected that shattering of 
rain droplets is enhanced through the combination of 
wet surfaces and larger droplets while shattering of pes¬ 
ticide drops is enhanced through lowering of surface 
tension. The bounce and shatter criteria developed in 
this paper assume a dry surface and additional work 
is required to characterise the influence of a wet sur¬ 
face on droplet impaction outcomes. The post-impact 
behaviour of droplets for both rain dispersal and pesti¬ 
cides needs to be characterised through further research 
to enable distribution of either pathogens or pesticide 
through a plant canopy to be determined. 


6 Conclusions 

The application of agricultural chemicals via spraying, 
and spreading of pathogens by rain splash, are compli¬ 
cated processes and there are many inherent variations. 
These include: the plant shape and structure with re¬ 
sulting variation in leaf angles; leaf microstructure that 
can contain, for example, hairs, trichomes, stomata, 
veins, and wax structures; variable droplet properties 
(in particular size, velocity, trajectories, density); vari¬ 
able formulation properties (e.g. surface tension and 
viscosity); variable droplet impact conditions (due to 
the increased possibility of impaction on a previously 
wet area as the run progresses); and variable airflow due 
to micrometeorology and turbulence resulting from the 
spray application. The experimental technique used in 
this study incorporated many of these variables, such 
as the use of real leaves with all their inherent natural 
variability, and commercial agricultural spray nozzles 
resulting in a range of droplet velocities, sizes and tra¬ 
jectories being produced. While there is scatter in the 
experimental measurements largely due to these uncon¬ 
trolled factors, the results do show broad agreement 
with predictions. In particular: 

— shatter of droplets increased with increasing size 
and velocity of droplets at impact; 

— shatter of droplets increased as the surface become 
harder to wet (shatter was more obvious on the hard 
to wet wheat leaves than the easy to wet cotton 
leaves); 

— for easy to wet surfaces or superspreader mixtures, 
the droplets mostly either adhered or shattered on 
impact, with very few droplets bouncing; 
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— bouncing of droplets was most pronounced on hard 
to wet surfaces (such as wheat) with high surface 
tension spray mixtures (such as water). 

Droplets which had a diameter of 200 /im or smaller and 
were travelling at less than 1-2 m/s mostly adhered on 
impact, with a few exceptions, particularly in the case 
of water droplets on wheat. On wheat, many of these 
smaller, slow moving droplets were observed to bounce. 
Shatter became predominant as the velocity and size of 
droplets was increased, without exception. Addition of 
adjuvants (surfactants) which lower surface tension can 
assist with initial adhesion, but will also lead to more 
shatter. By combining modeling with experimental ob¬ 
servations, a better insight into the processes occurring 
with different spray conditions has been achieved. 
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